Phase transition in a chain of quantum vortices 



00 
ON 
ON 

Or 
<D 



o 
o 

Oh' 

•4-3 

c3 



C. Bruder^B, L. I. Glazman 2,3 , A. I. Larkin 2,4 , J. E. Mooij 3 , and A. van Oudenaarden 3 
1 Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 
2 Theoretical Physics Institute, University of Minnesota, Minneapolis MN 554-55 
3 Delft University of Technology, 2600 GA Delft, The Netherlands 
4 L. D. Landau Institute for Theoretical Physics, Moscow 117940, Russia 

(February 7, 2008) 

We consider interacting vortices in a quasi-one-dimensional array of Josephson junctions with 
small capacitance. If the charging energy of a junction is of the order of the Josephson energy, the 
fluctuations of the superconducting order parameter in the system are considerable, and the vor- 
tices behave as quantum particles. Their density may be tuned by an external magnetic field, and 
therefore one can control the commensurability of the one-dimensional vortex lattice with the lattice 
of Josephson junctions. We show that the interplay between the quantum nature of a vortex, and 
the long-range interaction between the vortices leads to the existence of a specific commensurate- 
incommensurate transition in a one-dimensional vortex lattice. In the commensurate phase an 
elementary excitation is a soliton, with energy separated from the ground state by a finite gap. 
This gap vanishes in the incommensurate phase. Each soliton carries a fraction of a flux quantum; 
the propagation of solitons leads to a finite resistance of the array. We find the dependence of the 
resistance activation energy on the magnetic field and parameters of the Josephson array. This 
energy consists of the above-mentioned gap, and also of a boundary pinning term, which is differ- 
ent in the commensurate and incommensurate phases. The developed theory allows us to explain 
quantitatively the available experimental data. 
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I. INTRODUCTION 

The interest in Josephson junction arrays in the last 
decade was to a large degree prompted by the fact that 
these systems are suitable as a testing ground for var- 
ious predictions of quantum many-body theory (for an 
overview see, e.g., Refs. |l],|j). If the charging energy of 
a junction Eq = e 2 /2C is comparable with its Joseph- 
son energy Ej, the phase of the superconducting order 
parameter is subject to quantum fluctuations (here C is 
the capacitance of a junction). At some critical value of 
Ec/Ej the global phase coherence is destroyed, and the 
array becomes an insulatortj. This transition apparently 
is driven by proliferation of spontaneously created vor- 
tices, i.e., topological excitations of the array, in which 
the phase of the order parameter varies by 2tt on go- 
ing around a plaquette. At smaller ratios Ec/Ej vor- 
tices induced by an external magnetic field still possess 
quantum properties. The vortex dynamics is particularly 
sensitive to the quantum fluctuations of the phase: the 
vortex mass-pfor example, is finite entirely due to these 
fluctuations^. 

A single vortex in a Josephson junction array-,behaves 
as a ballistically propagating quantum particles. These 
particles are strongly interacting, however: for the val- 
ues of Ec /Ej < 1 at which the global phase coherence is 
preserved, their interaction energy U v - V oc Ej is larger 
than the bandwidth for a single vortex. A finite mag- 
netic field applied perpendicularly to the array, creates 
a lattice with a vortex density proportional to the field 
strength. Depending on the magnetic field flux per pla- 



quette, the vortex lattice is commensurate or incommen- 
surate with the junction arrayu. The commensurability 
effect exists of course even for classical vortices in an 
array with Ec/Ej — > 0. The array acts like a periodic 
potential with an amplitude U p ~ 0.2£j and some period 
a (the period of the Josephson array) for each vortexEx 
In a classical system, this is expected to be a source of 
strong pinning, as U p and U v - V are of the same order. 
Quantum fluctuations bring new physics into the prob- 
lem. The period of the pinning potential is relatively 
small, and therefore its amplitude is suppressed readily 
by quantum fluctuations. On the contrary, depending on 
the magnetic field, the vortex lattice period may be sig- 
nificantly larger than a, thus making the vortex lattice 
robust against quantum fluctuations. 

Commensurability effects in a chain of quantum vor- 
tices were investigated in the recent experiments of Oude- 
naarden et aln. There a number of two-dimensional ar- 
rays with various ratios Ec/Ej, and various widths of 
the order of 10 cells were studied. All arrays were quasi- 
one-dimensional in fact, their length varying between 100 
and 1000 cells. Superconducting contacts parallel to the 
long sides of the array were providing a potential confin- 
ing the vortices to the central row of the array, see Fig. [l]. 
Applying a current through the leads (perpendicular to 
the one-dimensional vortex chain) and measuring the re- 
sistance of the system as a function of the magnetic field, 
van Oudenaarden et al. found almost zero resistance in 
the regions centered around the commensurate values of 
the one-dimensional vortex density. This was interpreted 
as an indication of a finite-gap state ("Mott insulator" 
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phase), induced by the vortex- vortex interaction in the 
presence of a periodic potential. If the magnetic field was 
tuned away from these special regions, a transition to a 
resistive state was observed, indicating moving vortices 
( "conducting" phase) . 
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FIG. 1. Quasi-one-dimensional Josephson array. Each side 
of a plaquette corresponds to a single Josephson junction. 
Crosses denote vortices located in the central row of the ar- 
ray. The properties of the vortex chain are probed by passing 
a current / from one superconducting contact to the other 
and measuring the voltage between them. The array width is 
W, and the size of a single plaquette is a x a. 

The analogy between the observed transition in the 
system of vortices and the textbook Mott transition in 
electron systems can be made more explicit: Ec/Ej here 
plays the role of the ratio t/U of the electron bandwidth 
to the on-site repulsion potential; the magnetic field B, 
creating vortices, plays the role of the electron chemical 
potential fj,. The phase diagram in the variables (t/U, ji) 
consists of two phases. The insulating phase occurs at 
relatively small values of t/U < 1. In this phase, the 
electron density is constant as a function of fj,, and fixed 
by the commensurability condition (one electron per lat- 
tice site). One can assign the value fi = to the line 
of particle-hole symmetry in this phase diagram. Devi- 
ation from this symmetry line makes the excitation gap 
in the insulating phase smaller; the gap width is a non- 
analytical function of fi, reflecting the violation of the 
particle-hole symmetry at /i ^ 0. At a certain criti- 
cal value of which depends on t/U, the gap disap- 
pears and a transition to a "conducting" phase occurs. 
This conventional picture is modified somewhat in the 
one-dimensional case, but qualitatively remains valid. A 
similar description applies to Mott transitions in Bose 
systems with repulsiontj. 

The line of particle-hole symmetry in the case of vor- 
tices corresponds to a special value Bo of the magnetic 
field, which induces a vortex lattice commensurate with 
the period of the junction array. The analogy to the elec- 



tronic case described in the last paragraph suggests that 
the vortex density remains constant in a finite interval of 
B around Bq. One also expects a gap in the excitation 
spectrum of the array that diminishes as a function of 
\B — Bq\ within the commensurate phase (which is the 
analog of the insulating phase). Indeed, in the experi- 
ment a cusp-like dependence of the resistance activation 
energy on \ B— B \ was observed (see Ref. [|and Section[v| 
of this paper) . At some critical value of \B — Bq\ the ob- 
served temperature dependence of the resistance becomes 
considerably weaker. This may indicate that a transition 
to a gapless ("conducting") phase occurs, accompanied 
by the creation of discommensuration solitons. 

In this paper, we present a quantitative theory of the 
commensurate-incommensurate transition for a chain of 
quantum vortices in a quasi-one-dimensional Josephson 
array (see Ref. |ll| for another type of commensurate- 
incommensurate transition in Josephson arrays). The 
transition to the incommensurate state occurs by pro- 
liferation of the discommensuration solitons through the 
vortex chain. We are able to develop a theory by analyt- 
ical means, because of a special feature of the system we 
consider. It turns out, that the range of the interaction 
between the vortices is much longer than the inter- vortex 
distance. Therefore, the solitons consist of many vor- 
tices, and possess a large effective mass. Thus the theory 
for the commensurate-incommensurate transition is es- 
sentially classical. However, to relate the parameters of 
this theory to the generic properties (Ec and Ej) of the 
Josephson array, we need to consider a single vortex as a 
quantum particle: the amplitude of the periodic pinning 
potential depends on the bandwidth of the vortex. After 
that, we are able to find analytically the excitation gap 
existing in the commensurate phase and the boundaries 
of this phase in the (B ,Ec / Ej) plane. 

We relate the characteristics of the commensurate and 
incommensurate phases to an observable quantity, the 
activation energy of the resistance Er. In the commen- 
surate phase, the transfer of one flux quantum between 
the edges of the array occurs via a sequence of solitons 
propagating through it. The number of solitons neces- 
sary to transfer one vortex is equal to the ratio of the 
periods of the vortex lattice and the junction array; typ- 
ically this ratio is large. At any time during the vortex 
transfer, there is no more than one soliton present in 
the chain. We demonstrate that En depends not only 
on the properties of the "bulk" one-dimensional system, 
but also reflects boundary pinning effects, accompany- 
ing the passage of vortices through the ends of the array. 
One soliton changes the length of the vortex chain only 
by one period of the junction array, which is less than 
the inter-vortex spacing. Hence, in the commensurate 
phase, the process of vortex flow through the array can 
be viewed as motion of a rigid vortex chain. Because of 
the rigidity, the vortex chain cannot adjust itself to the 
boundary pinning potential. The potentials produced by 
the two ends of the array add to En: the relative phase 
of these two contributions depends on whether the total 
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flux piercing the junction array equals an integer number 
of flux quanta. Thus, in the commensurate state, there 
are two major terms in Er. The first term is the activa- 
tion energy of a soliton, and the second term is the sum 
of the boundary pinning energies. This second (smaller) 
term oscillates with the magnetic flux piercing the array. 
In the incommensurate state, the vortex chain is com- 
pressible, and can adjust to the boundaries of the array, 
if the latter is sufficiently long. As a result, the main 
term in Er is the boundary pinning potential, which a 
vortex has to overcome to enter the array. A correction 
is provided by the finite compression energy of the chain. 
Its average value depends on the compressibility of the 
chain, renormalized by the solitons, and is inversely pro- 
portional to the length of the array. This term oscillates 
with the flux threading the system. 

The paper is organized as follows. In Section 0, we 
introduce a model of classical vortices in the quasi-one- 
dimensional Josephson array. Here we neglect the dis- 
creteness of the array, and the screening of the vortex- 
vortex interaction. This approximation means that the 
the vortex chain is entirely incompressible. We establish 
the stability criterion for a one-dimensional vortex chain 
against formation of a zigzag structure. We calculate the 
equilibrium number of vortices as a function of the mag- 
netic field and determine the boundary pinning caused by 
the interaction of the vortices with the ends of the array. 
For an incompressible chain, this gives us the equilibrium 
position for each vortex. In the following part, Section 
[II , we discuss bulk pinning by reintroducing the discrete- 
ness of the junction array. The array creates a periodic 
potential for each vortex, which behaves as a quantum 
particle in this potential. We demonstrate that typically 
the amplitude of the quantum fluctuations of a vortex 
exceeds the period of the array. (Thisjustifies, in fact, 
the approximations made in Section ||). We calculate 
the residual pinning potential, suppressed by quantum 
and thermal fluctuations, acting on a sin gle vortex. 

The results of the Sections || and III are directly ap- 
plicable to short arrays, i.e., arrays that are shorter than 
the range of the vortex- vortex interaction. The main goal 
of these sections though is to provide us with the coeffi- 
cients necessary to write down the effective Hamiltonian 
describing a compressible chain in a long array. We start 
the next Section IV with an estimate of the vortex- vortex 



interaction range. It is defined by two mechanisms: (1) 
the effect of the magnetic field induced by the vortices, 
and (2) the interaction of the vortices in the Josephson 
array with the Abrikosov lattice in the contacts to the ar- 
ray. The estimate|-demonstrates that for the conditions 
of the experimentsa the range indeed exceeds greatly the 
inter-vortex distance, but still may be smaller than the 
system length, making it necessary to account for a finite 
compressibility of the vortex chain. We therefore derive 
the long-wavelength theory for the compressible vortex 
chain. This theory enables us to describe, in Section [v|, 
the commensurate-incommensurate transition. We de- 
termine the boundaries of the commensurate phase, and 



find the dependence of the activation energy for elemen- 
tary excitations on the parameters of the system. Also in 
this section, we discuss the behavior of the resistivity fol- 
lowing from the picture we developed. We compare our 
results with the existing experiment in Section 
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Us- 
ing the experimental values of the vortex density at the 
commensurate-incommensurate transition and the max- 
imum of the activation energy of the resistance, we are 
able to give parameter-free estimates of the range of the 
vortex-vortex interaction and of the elastic constant of 
the vortex chain. The effective pinning potential turns 
out to be at least an order of magnitude smaller than the 
bare potential due to quantum fluctuations (as calculated 
in Section |l|). The soliton length is extremely large and 
of the order of the length of the array. The long-range 
nature of the vortex-vortex interaction leads to a large 
value of the elastic constant: the chain is virtually rigid 
in the incommensurate phase. Our theory explains con- 
sistently the main experimental observations reported in 
in this paper and in Ref. g (1) the cusp-like dependence 
of the activation energy on the magnetic field in the com- 
mensurate phase; (2) the large value of this activation 
energy (compared to Ej and Ec), and (3) oscillations 
of the resistance with the applied magnetic field, with a 
period corresponding to one flux quantum through the 
entire array, in the incomme nsur ate phase. We conclude 
with a discussion in Section VII. 



II. RIGID VORTEX CHAIN 

We consider a two-dimensional Josephson array of lat- 
tice constant a, length L and width W where L ^> W, 
see Fig. |]. The "sites" of this array are superconducting 
islands, linked by Josephson junctions that are character- 
ized by a capacitance C and a critical current Ic- The 
phases ip\ of the order parameter of the islands (num- 
bered by vectors i) are the only dynamical degrees of 
freedom of the system. For an infinite two-dimensional 
system, the Lagrangian can be written in the standards 
way, 



£ 

<i,j> 



h 2 ( d<pi\ 



8E C V dt 



-Ej [1 - cos (pij)] 



(1) 



Here the sum is taken over the nearest neighbors, and 
c/Jij is the phase difference across a link of the array, 
Ej = Ic&o/Zk an d Ec = e 2 /2C are the Josephson 
and charging energy, respectively; $o = hc/2e is the 
flux quantum. The Lagrangian (Q) describes quantum 
fluctuations of the phase in the array. At a certain crit- 
ical value£3 of the ratio Ec / Ej ~ 1 , the proliferation of 
spontaneous vortices and antivortices through the system 
destroys the long-range order. We consider smaller val- 
ues of Ec/Ej, and neglect the existence of spontaneous 
topological excitations. Vortices in the spatial distribu- 
tion of the phase ip are then induced only by an external 
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magnetic field B. A vortex is characterized by a phase 
change of 2n on going around a plaquette. The effective 
Lagrangian in terms of the vortex positions, 



M 



dt 



(2) 



can be derived! from Eq. (|l|). Here M = ir 2 ti z /4a 2 E c 
is the vortex mass, U v — v (ri,Tj) is the interaction energy 
between the vortices, and U p (ji) is the pinning poten- 
tial which represents the effect of a discrete lattice of 
junctions on the vortex motion. In an infinite array, the 
energy U v - V depends only on the distance between vor- 
tices, and can be approximated by the standard expres- 
sions-valid for vortices induced in a thin superconducting 
filmEx For a geometrically restricted array, Fig. [|, the 
form of the interaction potential U v - V depends crucially 
on the boundary conditions for the phase that are set by 
the massive superconducting contacts. The superfluid 
density in these superconducting strips exceeds greatly 
the effective superfluid density in the array. Therefore, 
each vortex in the array is repelled from the boundaries 
(this is represented by the terms J7„— u( r i) r i) = Uv-v{yi) 
in the Lagrangian). At a sufficiently weak magnetic field, 
B &0/W 2 , the inter- vortex distance is large enough, 
and vortices occupy only the central row of the array. 

The pinning potential in this one-dimensional case may 
be modeledcl by a function of x only, 



U p (x) = 0.1Ej[l - cos(2irx/a)} 



(3) 



We will see in the next section, that quantum fluctuations 
of the vortex positions strongly diminish the role of pin- 
ning by the periodic lattice of Josephson junctions. For 
now, we will ignore the contribution of the pinning poten- 
tial, given by the third term in the Lagrangian Eq. (Q) . 

Because of the large superfluid density in the contacts, 
the phase of the order parameter varies only slightly 
along each of the long boundaries of the array. Currents 
induced by a vortex in the array flow through it almost 
perpendicularly to the boundaries. In the limit of infinite 
superfLuid density and infinite London-Pearl penetration 
depthEj in the contacts, the currents within the array do 
not decay with the distance from a vortex. As a result, 
the range of the vortex- vortex interaction is infinite. The 
interaction potential, up to an arbitrary constant, has the 
following form 



-2w 2 Ej\ 



,\/W 



(4) 



at | afj ~ Xj \ ft W. Because of this form, the vortex chain 
is absolutely rigid at small wave vectors. 

If the smallest inter-vortex distance exceeds the array 
width, one can use the limiting form of the potential (^) 
to calculate the contribution of the vortex-vortex inter- 
action [the second term in Eq. (^)] to the energy of the 
vortex chain. It is more convenient, however, to write 



down this energy directly in terms of the phase distribu- 
tion in the array: 



dtp 1 



N 



tt sign(x — Xi) + ifo] 



(5) 



Here, we have replaced the phases ip\ of the islands by 
a continuous variable <p(r). The form (||) of the phase 
gradient is valid at distances | ac — rc» ] > W away from 
the vortex centers x,. The magnetic field enters via the 
one-dimensional density n = BW/ $0 ■ The phase (po has 
the meaning of the average phase difference between the 
contacts, and will be used as a Lagrange multiplier to 
enforce the condition of fixed current I in the y-direction 
through the array. The typical shape of the phase gradi- 
ent is illustrated in Fig. |^. The energy in the presence of 
a current / between the contacts can be written as 



E({xi}, ipo) 



EjW 



-l/2 \oy 



f>. (6) 




FIG. 2. Current (in arbitrary units) across the array as a 
function of the coordinate x along the array. Note the jumps 
at the vortex positions Xi. For nL = integer (solid line), we 
can identify the two edges at ±L/2, and shifting the vor- 
tex chain along the x-direction does not change the energy, 
Eq. (^). In contrast to that, for nL 7^ integer (dashed line), 
the energy depends on the position of the chain (boundary 
pinning) . 



The equilibrium positions of the vortices and tpo for 
a given value of the current I are defined by the set of 
conditions 



dE({xi},ip ) 

dxi 
dE({xj},ip ) 

dipo 



= 
= . 



(J) 



At first we will consider the case 1 = and will de- 
termine the equilibrium number N and positions X®, 
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i = of vortices in the array. Solving Eqs. (0), 

we obtairO 



N = Int(nL) 



(8) 



where Int(s) is the integer part of x. We will consider 
only positive values of the magnetic field, n > 0. The 
equilibrium positions are given by 



2i - 1 - N 
2n 



(9) 



which means that the vortices are equidistant, x\ 
x\ = 1/n. The first and the last vortex of the chain are 
located at a distance [nL — Int(ni) + l]/(2n) > l/(2n) 
away from the ends of the array. On increasing the flux, 
they move towards the center. 

For deviations of the vortex coordinates from their 
equilibrium positions, the energy Eq. (^J) may be ex- 
pressed as 



E — En 



2Tr 2 Ejn 



W 



N 

{5 2 Frac(?ii) + ^(xj - x°) 2 } , (10) 



where x = Xi /N is the center-of-mass of the vortex 
chain; for brevity, hereafter we use the notation 



Frac(ni) = nL — Int(nL) 



(11) 



Each individual vortex resides in a parabolic well, cen- 
tered at the vortex equilibrium position; this is the result 
of the infinite-range interaction between the vortices. 

The term proportional to x 2 is caused by the interac- 
tion of the vortices with the two boundaries at ±L/2: 
if nL is integer, shifting the vortex chain along the x- 
direction does not change the energy, Eq. (||), see Fig. ||. 
For general values of nL, the energy depends on the po- 
sition of the chain. That means that the boundaries pin 
the vortex chain. 

The activation energy of the system is given by the 
difference in ground-state energies Eq of the N + 1 and 
N- vortex chains at a given value of the flux density n. A 
straightforward calculation starting with Eq. (^|) at / = 
yields 



E h (nL) = ^-^-[1 
V ; 2Wn l 



Frac(ni)]Frac(nL) 



(12) 



The boundary pinning energy (|12|) vanishes for inte- 
ger values of nL; the maxima between two zeroes are 
w 2 Ej/(8nW), and decay with the magnetic field as 1/n, 
see Fig. H 




FIG. 3. Activation energy (in units of tt 2 E.jL/2W) of the 
rigid chain as a function of nL = BWLf&o. The critical cur- 
rent is given by the same curve: J^ ray — 4_Eb(nL)/$o, see 
Eq. @ and also Ref. jT|. 



It is also possible to calculate the critical current of 
the array, L^ ray . We define it as the current at which 
the stability of the center-of-mass x is lost. The relation 
between J and x can be found from Eq. (^j, 



I = 



2irE 



$o J-L/2 \oy 
4ir 2 Ej x 



Frac(nL) 



(13) 



The center-of-mass stability requires that no vortex is to 
enter or leave the system; this restriction leads to 



rarray_ 
1 C 



4£ b 
$0 



nW 



1 - Frac(nL)]Frac(nL) . (14) 



This functional dependence of the critical current on nL 
for the array coincides with the one obtained in Ref. [l5| 
for a thin-film bridge. 

It is clear that the one-dimensional approximation 
breaks down for large magnetic fields, i.e., if the vortex- 
vortex distance becomes considerably less than the width 
W of the array. To obtain a quantitative value of the crit- 
ical field, we studied the instability of the vortex chain 
towards formation of a zigzag deformation. For that pur- 
pose, we derive the full formula for U v — v (ri,Tj), which 
requires properly taking into account the infinite number 
of image vortices necessary to fulfill the boundary condi- 
tions at the superconducting contacts. After that, we re- 
place Eq. (0) by its two-dimensional analogue, which de- 
pends on the two-dimensional vectors of displacements of 
each vortex. An analysis of the dependence of this energy 
on the transverse vortex displacements yields the value 
of the critical field at which the zigzag pattern forms: 



B CT it — 



0.65W 2 



(15) 
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In other words, the vortex distance has to be larger 
than 0.65VF for the one-dimensional approximation to 
be valid. 



III. PINNING BY THE PERIODIC POTENTIAL 

In the last section, the discreteness of the system 
was neglected completely. At first sight, this seems 
to be an unreasonable approximation. Indeed, only at 
nW > (0A/ir 2 )(W/a) 2 the amplitude Q.2Ej of the pin- 
ning potential (||) is smaller than the variation SE of the 
energy ( |Io|) if a single vortex is displaced by a/2 (a/2 is 
the distance between the minimum and maximum of the 
pinning potential) . Note that for the stability of a single- 
row vortex chain, the condition nW < 1.62 must be sat- 
isfied (see Eq. (|lq)). The two restrictions on nW are 
incompatible except for quite narrow arrays, W/a < 6. 
In this section, we demonstrate that the effective pinning 
potential is reduced significantly by quantum fluctuations 
of the vortex coordinates, which makes the above restric- 
tion irrelevant, even at relatively small ratios Ec/Ej. 

Forgetting the interactions with the other vortices for 
a moment, each vortex is described by a Hamiltonian 



H= 2M + 0AEj ^ ~ C0S ( 2 W«)] 



(16) 



That means, it is a delocalized quantum particle char- 
acterized by a band structure e(k). In the limit 
of small quantum fluctuations of the phase, Ec ^$ 
OAEj, the problem can be treated in the tight-binding 
approximation!!! This leads to the following expression 
for the dispersion relation 



E 

e(k) = — ^- cos(fca) , 
where the bandwidth is given by 
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(17) 



E p = -y/0.1EjE c exp(-2y/0.1Ej/Ec) , (18) 

and the effective mass of the vortices in the periodic po- 
tential by 



a 2 E n 



m 



eff 



(19) 



At stronger fluctuations the tight-binding approximation 
is inadequate, and the bandwidth becomes of the order of 
Ec- We note that E p can be interpreted as a transition 
amplitude. The exponent of this amplitude can be also 
extracted from Ref. ||, where the rate of transitions be- 
tween two adjacent minima of the pinning potential was 
estimated. Ac cording to Ref. [|, this exponent is approx- 
imately 2.25a/ Ej/8Ec which is about 10% higher than 
the exponent in Eq. (fl8|). 

In Eq. ( ]To| ) we showed that each vortex moves in a 
parabolic potential produced by the interaction with the 



other vortices. The corresponding oscillation frequency 
of a particle having effective mass m e ff is 



Uir 2 Ejn _ 1 

Wm e ff h 



2-rT 2 EjE p na 2 /W , (20) 



and the mean-square oscillation amplitude is 



((Xi 



WeffWosc 



E p W 

8ir 2 Ejn 



(21) 



The quantum fluctuations implied by Eq. ( pl| ) lead to 
a reduction of the effective pinning potential, which be- 
comes now a (periodic) function of x®. Estimating the 
effective pinning, we assume that the inter-vortex inter- 
action on the scale of the lattice constant a is weak com- 
pared to E p . In accordance, with the standard prescrip- 
tion of solid-state physicscj, we replace the quasi-wave 
vector in the dispersion relation Eq. ( p"7| ) by an operator 
p/h, and consider the Schrodinger equation for a quasi- 
particle with Hamiltonian 



H 



2 / - 0\ 2 

-m off w osc (a; - xj 



(22) 



Here p and x are canonically conjugate variables. But 
now we can view x as the momentum of some particle, 
moving in the potential e(p/h) that is periodic in the 
coordinate p of the particle. Therefore, eigenstates of 
the Hamiltonian ( [22] ) with various values of the "quasi- 
momentum" X® form bands. For each band, the energy 
is a periodic function of x® with period a. At zero tem- 
perature, we are interested in the lowest band, with the 
energy 



Uf{x1) 



Uf cos(27ra£/a) 



(23) 



The value of Up S depends on the magnitude and form of 
the periodic potential e(p/K). Using (17) and (|i~8|), we 
find: 



Uf = Ej ^2nW§-^ 



2V2W [EL 1 



Ej nW 



(24) 

Equation (^3|) gives the effective pinning potential for 
a single vortex. It is worth noting that the pinning 
strength diminishes with the increase of the equilibrium 
inter- vortex distance l/n. 

The approximations we employed in deriving the form 
[Eq. (p3[)] and amplitude [Eq. (pi]) ] of the pinning po- 
tential require a sufficiently wide band for the motion of 
a vortex in the periodic potential. In other words, the 
exponential factor in Eq. must be small. In the op- 
posite limit of negligible quantum fluctuations, the mag- 
nitude of the effective potential is 0.2Ej, and the function 
Up S (%i) has cusps at x° coinciding with the maxima of 
the bare potential U p (x) defined in Eq. (||). Each cusp 
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in C/p ff (x") corresponds to a jump of the coordinate of 
a classical vortex between the minima of the potential 
U p (x). 

One may get an idea of how effective the quantum 
smearing is, by estimating Up S at W/a = 10 and 
Ec — OAEj which is close to the limit of applicabil- 
ity of ( |l8| ) . Substitution of these values in Eq. (g4j) yields 
Uf ~ 0m6y/nWEjexp(-3.58/VnW). In this exam- 
ple, the effective pinning potential gets smaller than its 
bare value at 1/n > W/40 which is always the case in 
practice. 

At finite temperature T > huj osc , we have to consider 
the averaging of the periodic potential by quantum and 
thermal fluctuations; these further reduce the pinning. 
Summing the geometric series, we obtain 



uf(T) = u; 



1 + exp(-huj osc /T) 



(25) 



We will now calculate the pinning of the rigid vortex 
chain. Each vortex is subject to the potential Eq. (23). 



Summing over the members of the chain [which are lo- 
cated at the positions shifted by x from the equilibrium 
values CI9|)1 leads to 



/v 



U pin (x) = J2uf(x°i+x) 




0.01 



n-n„ 



FIG. 4. Activation energy entering the resistance as a func- 
tion of n = BW/&o- The energy plotted here is the sum of 
the boundary pinning term Eq. (^2|), and twice the absolute 
value of the bulk pinning term Eq. (|26[). Energy units as in 
Fig. ^, and n, no are measured in units of 1/a. 



Although each vortex is a quantum particle (as we have 
stressed at the beginning of this section), the possibility 
of vortex permutations may be safely neglected: for each 
vortex, ((xi — x°) 2 ) <C n~ 2 , i.e., the oscillation amplitude 
is much less than the inter-vortex distance. 



-Uf{T) cos(27ra/a 



sin(iV7r/na) 
sin(7r/na) 



(26) 



For commensurate values of the flux, i.e., if l/na is 
integer, we get 



U pin (x) = -Uf(T)cm{2nx/a)N ; 



(27) 



the pinning barrier is proportional to the total number 
of vortices, i.e., the pinning is strong. In the immedi- 
ate neighborhood of the commensurate points, however, 
there are values of n for which 



Int(nL) 



(28) 



At these vortex densities, the numerator of Eq. (|26|) van- 
ishes, i.e., there is no pinning. The spacing between these 
zeroes is approximately given by na/ L, which may be less 
than 1/L. The rapid oscillations are caused by the fact 
that we are considering a completely rigid vortex chain. 
If we neglect the oscillations, and just look at the maxima 
of J7pi n , it turns out that the pinning strength behaves as 
l/\n— 7io| close to commensurate densities Uq. 

The activation energy for the resistance can be esti- 
mated as the sum of the amplitude of U p i n (x), Eq. (26), 
and the boundary pinning term E\>(nL), see Eq. (12). 
The result is shown in Fig. |i[ 



IV. COMPRESSIBLE VORTEX CHAIN 

In the last section we considered the case of the 
infinitely long-range vortex-vortex interaction with the 
consequence that the vortex chain was completely rigid. 
We will now discuss the importance of screening and the 
resulting compressibility of the chain. 

Screening of the vortex-vortex interaction in the 
Josephson array is due to two effects: (1) screening by 
the magnetic field created by currents flowing around the 
Josephson vortices (Meissner effect), and (2) interaction 
with the vortex medium in the contact pads. To start 
with, we consider the first of these two effects. 

The distribution of currents flowing around a vortex 
depends on the dimensionality of the system. The cur- 
rents around a vortex line in a three-dimensional super- 
conductor drop off exponentially fast, the characteris- 
tic length being the London penetration depth Al- For 
a vortex in a superconducting, film of small thickness 
s <C Al, the screening lengthEZI is thickness-dependent, 
A = A|/s. At a sufficiently large distance from the cen- 
ter of a single vortex, r > A, the spatial distribution of 
these currents is controlled by the Meissner effect, 



47r 2 r 2 



(29) 



Note, that the distribution ( p9| ) is insensitive to the short- 
scale structure of the two-dimensional vortex: it may be 
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a Josephson vortex in the array, as well as an Abrikosov 
vortex in the contacts. In both cases, the currents in- 
duced by a vortex fall off as 1/r 2 at sufficiently large 
distances r from its center, resulting in a vortex-vortex 
interaction with a finite range A s . In the specific case of 
a quasi-one-dimensional Josephson junction array con- 
tacted by superconducting films, we may estimate A s by 
matching the current density j\ with the density ji of 
the current flowing around a vortex in the array, 



J2 



2it 2 Ejc 
<J> W 



(30) 



The current density ji corresponds to a single- vortex con- 
tribution to the phase gradient (^). Equating ji (A s ) w j% 
leads us to the estimate 



As 
W 



8tt 4 EtW 



(31) 



At distances r > A s , the anisotropy of the system is not 
important for the current distribution, and Eq. (|2^) is 
applicable. At smaller distances, the current distribu- 
tion is highly anisotropic. In the cc-direction, the screen- 
ing length is A s , while in the y-direction the currents are 
confined to the thin-film penetration depth A, which de- 
pends on the properties of the leads. Typically the sheet 
superfluid density p s in the contacts exceeds greatly the 
effective superfluid density pj in the Josephson junction 
array, which leads to A -C A s . 

The electrodynamic effect limiting the radius of inter- 
action considered above exists for any geometry of the 
contacts. The estimate (|3l]) does not depend on the char- 
acteristics of the contact material. However, the numer- 
ical coefficient in (|3l]) is geometry-dependent. In Section 
VI we will further refine tha^estimate ([H]) for the specific 
geometry of the experiment!]. 

In the above consideration, we have completely ignored 
the existence of Abrikosov vortices in the contacts. This 
is acceptable only if their currents do not overlap with 
the currents created by Josephson vortices. The two cur- 
rent density fields are spatially separated if the distance 
d of the last row of the Abrikosov lattice to the edge of 
the junction array (see Fig. [j]) exceeds A. In the opposite 
case d < A, the Abrikosov lattice effectively truncates 
the currents created in the contact by a Josephson vor- 
tex. Indeed, a small shift of the lattice in the direction 
perpendicular to the edge of the array is sufficient to com- 
pensate these weak currents. To estimate the interaction 
potential range A s in this case, it is sufficient to deal with 
the energy of the supercurrents, 



E 



4:ir 2 Ej\ s 
W 



16A S A 



(32) 



and neglect the magnetic fields the supercurrents create. 
In Eq. (|3^), the first term corresponds to the energy of 
currents in the array that flow in the region \x\ < A s 



around the vortex. The second term is the energy of the 
supercurrents in the contacts. These currents are trun- 
cated at the position d ~ ^JWq/B of the first cow of the 
Abrikosov lattice, which numerically turns outEa to be a 
sound approximation. 

Minimization of the energy Eq. (|32|) with regard to A s 
yields: 



As 
W 



32ir 2 EjW A 



(33) 



This expression is valid for d < A, and at d ~ A it rea- 
sonably well matches the estimate (|3l| ). For typical ex- 
perimental values, Eqs. ( |3l| ) and (|33|) yield a screening 
length for which rT 1 <C A s < L. In other words, the 
vortex-vortex interaction has long, but finite range, and 
the vortex chain is not completely rigid. 

We will now develop a continuum description of the 
compressible vortex chain, i.e., we will express the energy 
of the chain in terms of the deviations u(x°) = x t — x®, 
and then go over to a deformation field u(x®) — > u(x). 
The energy of the chain will be the sum of a bulk pin- 
ning term and a boundary pinning term as before. In 
addition to that, there will be an elastic energy term. 

It is straightforward to express the bulk pinning term 
in terms of u(x): 

N 

E/pi„ = -[/ p off ^cos(27r^/a) 
i=l 

■ L/2 



f 1 2tt 
U% \ efacos[ — (u + ax)], 
J-L/2 a 



where 









a = na 


— - Int 






na ' 





no — n 
n 



(34) 



(35) 



is a dimensionless measure of the deviation from the com- 
mensurate value l/noa = Int(l/na). 

The ends of a long array (L> A s ) act on the compress- 
ible vortex chain as two independent sources of boundary 
pinning. In the case of an almost rigid chain, A s 3> 1/n, 
we can find the pinning potential created by a single end 
(say, the one corresponding to x — L/2) by a slight mod- 
ification of Eq. (0) . Namely, we introduce an exponential 
factor exp[7(x — L/2)] into the integrand, and replace 
the lower limit of integration by — oo. After that, we find 
the extremal value of ifo as a function of the position of, 
say, the last vortex in the chain xjv, and take the limit 
7 — * +0. This procedure yields 



U h {u) 



7T 2 E 



J 



[4(? 



(36) 
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Here for convenience we have introduced a new variable 
u instead of the coordinate xn, 

'L_l_ 
2 2n, 



xn 



(37) 
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and the coordinate x^ is no more than one half-period 
away from the end of the array, |u| < l/2n. The top of 
the boundary barrier is at nu = — 1/V12, and its ampli- 
tude is approximately 0.64Ej /nW. 

The elastic energy of vortices, which interact by long- 
range forces, in the long-wavelength limit takes the form 



L/2 dx (— 

-L/2 \9x 



(38) 



The elastic constant K can be expressed through the 
vortex-vortex interaction potential as 



K = n' 



dxU v - v (x) 



(39) 



Note, that the potential U v _ v (x) here is defined differ- 
ently from Eq. (g). Unlike in Eq. (|), we remove the 
uncertainty in the definition of the potential by requiring 
U v - V (x — > ±00) = 0. 

In order to estimate K, we adopt the following model 
for the interaction potential: 



U v - V [x)= A s exp 

W \ A s 



which correctly reproduces the cusp [cf. 
\x\ <C A s , and reaches zero at \x\ — > 00. 
model, we find 



(40) 



Eq. (§)] at 
Within this 



K 



dxU v - v {x) 



W 



-{n\ s f 



(41) 



(The real interaction potential in the planar geometry of 
the contacts considered above falls off as x~ 2 , rather than 
exponentially. However, this should not significantly al- 
ter the estimate) . We complete the estimate of K by 
adopting Eq. (31) for the value of A s , which yields 



K 



<I>Q?1 2 



(42) 



The elastic constant becomes softer, if the period of the 
Abrikosov lattice in the contacts is smaller than A, see 

Varying E = U e i + U p i n with respect to u(x) leads to 
the static sine- Gordon equation 



nUp S sin 



2nu 
a 



ax 



(43) 



This equation has been studied in many contexts, e.g., 
commetisurate-discommensurate transitions in adsorbate 
layerslia, here a is the difference of the lattice constants 
of the substrate and the adsorbate. Another example is 
the theory of long Josephson junctionstj where a is pro- 
portional to the magnetic field threading the junction. 



V. PHASES 

The behavior of a one-dimensional vactex chain is 
closely related to that of an adsorbate layerO: if the mag- 
netic field is commensurate, n = no, the vortex chain is 
commensurate with the junction array. The activation 
energy of an elementary excitation at n — no is given 
by the energy to push one soliton into the system. The 
length of such a soliton is given by 



K 

nUf 



and its energy is 



E'„ = — jKnUf 



(44) 



(45) 



A comparison of x s with the interaction radius Eq. ( |3l| ) 
yields x s /A s ~ (a/W)VnW J 0An 2 Ej/U£ s ; here we have 



{a/wy 

used the estimate ( p2[) for the elastic constant K. The 
applicability of Eq. (|43|) requires x s > A s , and therefore 
Eqs. (|44|), ( |45"| ) are valid only if the pinning potential is 
reduced by quantum fluctuations compared to its classi- 
cal value. 

On moving n away from no, the magnetic field tries 
to enforce a period of the vortex lattice that is different 
from the period of the pinning potential. The chain stays 
locked to a commensurate state up to a critical value of 
\n— no\, or, in other terms, until \a\ is less than some crit- 
ical value ac- Below the threshold, at \n — no| < no|ac|j 
the activation energy will diminish linearly with increas- 
ing \a\. This can be seen immediately from the analo- 
gous situation in a long homogeneous Josephson junction 
where quantized fluxons play the role of solitons. The en- 
ergy to create the first fluxon in the junction has some 
value at H = 0, and decreases linearly like —M\H\ with 
\H\; here M > is the magnetization of a single fluxon, 
i.e., a constant. In a complete analogy with this, in our 
case the energy to create a soliton in the commensurate 
phase is of the form 



K(n) =E S 



1 



no 



n a c 



(46) 



Note that the energy E s {n) has a cusp-like dependence 
on n — no. The point n = no is special: the creation 
of a soliton or anti-soliton costs the same energy. This 
situation is similar to the line of particle-hole symme- 
try for a Mott insulator. Deviation from the symmetry 
point makes creation of solitons, or anti-solitons prefer- 
able. This violation of the symmetry is the origin of the 
non-analytic dependence of E s on n. The critical values 
of a at which the soliton energy turns zero are given by 



ac = — 1 

7T 



InUf 
K 



(47) 
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Above the threshold, at \a\ > ac, discommensurations 
will exist: the chain is strained in the discommensura- 
tions, but this is offset by the fact that the rest of the 
chain can stay in the minima of the pinning potential. In 
this incommensurate phase, the concentration of solitons 
is finite. Due to the solitons, the vortex chain regains a 
finite compressibility K s , which depends on how far the 
system is tuned away from the critical points a = ±ac- 
Without giving the details here, we note that the depen- 
dence of the renormalized elastic constant on the control 
parameter a can be presented in a parametric formEj, as 
follows: 



Ks 
K 

H 

ac 



7T 2 £[l/ 7 #(7)] 
E(rt) 
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(48) 



Here if (7) and E{~{) are the complete elliptic integrals of 
the first and second kind, respectively. The chain softens 
near the critical points, where the proper expansiorO of 
Eq. (Ea) yields: 



Ks 
K 



8 Id 



a c 



ac 



hi- 



ac 



a 



ac 



(49) 



The softening occurs, because the solitons in the chain 
are rare, and the pair potential acting between them is 
exponentially small, U s ~ E s cxp(— x/x s ). Far away form 
the transition, at \n — no\ 3> ac^o, the solitons overlap, 
and K s = K . 

A finite voltage between the contacts to the array (see 
Fig. |l|) is related, by the Josephson relation, to the av- 
erage velocity of vortices moving along the array. The 
transport of a vortex through the system can be viewed 
as propagation of solitons through the vortex chain. The 
availability of solitons in the chain will clearly affect the 
resistance of the array. In the commensurate phase, the 
soliton density is exponentially small at low tempera- 
tures, E s (n) being the corresponding activation energy. 
In the incommensurate phase, there are mobile solitons in 
the system even at zero temperature. If one neglects the 
existence of boundary effects, the activation energy Er of 
the observable quantity, viz., resistance, would coincide 
with the activation energy of a single soliton. Hence, one 
would expect En = E s (n) in the commensurate phase, 
and En = in the incommensurate phase. 

It turns out, however, that boundary pinning modifies 
this picture. First of all, it may affect the ground state 
of the vortex chain. In the incommensurate phase, even 
weak boundary pinning will lead to a deformation of a 
long compressible chain. In the commensurate phase, the 
structure of the ground state starts to depend on the ra- 
tio of the boundary pinning energy Eh, ~ Ej/nW, and 
the soliton energy E s (n); this ratio depends on the bare 
parameters of the system, and may be small or large. 
Second, the set of excited states the chain goes through 



during the elementary act of a vortex transfer, also de- 
pends on the boundary pinning. These two factors deter- 
mine the dependence of En on the characteristic energies 
E s and E\>. We will analyze the activation energy En for 
both cases of small and large value of this ratio. 
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FIG. 5. Phase diagram of the compressible vortex chain: 
activation energy as a function of n — nc- (a) Eh S> E s , i.e., 
boundary pinning dominates (E s — 0.5). (b) Eh -C E s , i.e., 
soliton formation energy dominates (E s = 5). Energy units as 
in Fig. ^, and the one-dimensional vortex densities n, no, nc 
are measured in units of 1 /a. On the incommensurate side of 
the transition, n > nc, solitons will form spontaneously. The 
physics of the incommensurate phase is therefore determined 
by boundary pinning and the elastic energy and is the same 
for (a) and (b). 



In the commensurate phase, and in the presence of 
strong boundary pinning, E\, > E s (n)/2na, the vortex 
chain in the ground state will adjust itself to the length 
of the array to minimize the pinning energy. This means 
there are solitons in the ground state, unless nL is an in- 
teger. The largest number of solitons in the ground state 
occurs at a half-integer value of nL, and equals l/2na. 
In this particular case, the chain without solitons is the 
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configuration with the highest energy that the system 
passes through during a vortex transfer. In this state, 
the chain is not adjusted to the boundary of the array, 
and the boundary pinning energy attains its maximum 
value i?b- The difference of this energy from the ground 
state is Er = E^ — E s (n) /2na. In the opposite case of in- 
teger nL, there are l/2na solitons in the "saddle point" 
state, and the activation energy reaches its maximum, 
Er = E h + E s (n)/2na, see Fig. |a. 

If -Eb < E s (n), there are no solitons in the ground 
state of a commensurate chain. Moreover, during the 
process of a vortex transfer through the array, there is 
at most one soliton in the chain. Since a soliton changes 
the length of the vortex chain only by a, the chain re- 
mains rigid on the scale of the inter- vortex distance l/n. 
Therefore, we arrive at the following picture of the vortex 
transfer. The passage of each soliton shifts the chain by 
a. The transfer of a vortex requires the sequential pass- 
ing of 1/na solitons. In this process, the chain moves as a 
rigid object in the presence of boundary pinning. Thus, 
Er is the sum of E s (n) and the boundary pinning energy 
Eq. (O) for a rigid chain, 



E R = E s (n) + E h (nL) 



(50) 



see Fig. ^d. The soliton energy pq ) vanishes at the 
boundaries of the commensurate phase. Before it van- 
ishes, we cross over to the case described in the previous 
paragraph. 

We will now discuss the incommensurate phase, n > 
nc- At the phase transition, the soliton formation en- 
ergy vanishes, and solitons will start to form sponta- 
neously. Correspondingly, the physics of the incommen- 
surate phase will be determined by boundary pinning 
and by the elastic energy, and the behavior of the activa- 
tion energy will be identical for the two panels of Fig. ||. 
The chain is compressible; the elastic constant K s (n) is 
renormalized down by solitons, see Eq. ( f48| ). The ad- 
justment of the vortex chain to the length of the array 
leads to a finite deformation. The corresponding elastic 
energy can be found with the help of Eq. ( |3S| ) with K 
replaced by K s (n). The maximum value of the defor- 
mation du/dx = l/2nL corresponds to half-integer nL, 
and the elastic energy associated with it is K s (n)/8n 2 L. 
For large L, this energy is inevitably smaller than the 
boundary pinning potential (|3^). To initiate an elemen- 
tary act of vortex transport through the array, a shift 
of the end vortex through the maximum of the potential 
( |36| ) should occur. This varies the deformation of the 
chain by l/ny/3, or by (1 - l/Vfy/n. The corresponding 
elastic energy in both cases is the same, and is equal to 
K s (n)(l - 2/v / 3) 2 /8n 2 Z,. The net variation of the elastic 
energy involved in the described shift of the vortex, is 



5U e 



K s (n) 
8n 2 L 



1 



2 

71 



i 



-0.12 



K s {n) 
n 2 L 



(51) 



the boundary pinning amplitude. At some other values 
of nL, the variation in the elastic energy involved in the 
process of passing the boundary barrier, attains its maxi- 
mum value —5U C \. Thus, the resistance activation energy 
oscillates between two values, 



Er 



Ei KJn) 
0.64— J - ± 0.12-^/ 

nW n A L 



(52) 



with the period An = 1/L. 



To obtain the activation energy for the resistance Er 
at this particular value of nL, one should add 5U C \ to 



VI. COMPARISON WITH THE EXPERIMENT 

The resistance R of a number of arrays of Joseph- 
son junctions was measured in the geometry depicted in 
Fig. |l| in the presence of a magnetic field. Arrays with 
lengths L varying between 100a and 1000a, and widths 
W of 7a and 3a were studied. The characteristic Joseph- 
son energy for all the samples was about IK, with the 
ratio Ej/Ec varying within the limits 0.7 to 2.8. (The 
details of sample preparation as well as the experimen- 
tal techniques can be found in Ref. |). The main qual- 
itative feature of the field dependence of R consists in 
the existence of a finite region of magnetic flux densi- 
ties n around the commensurate value no ~ l/3a, where 
the resistance is strongly suppressed (Mott phase for the 
system of quantum vortices). The width of this region 
becomes smaller with the increase of the "quantum pa- 
rameter" Ec/Ej, see Ref. ^, in agreement with the no- 
tion of the Mott transition. 

Within the Mott phase, the resistance clearly displays 
an activated behavior, with the activation energy Er 
strongly depending on the deviation \n — no I from the 
point of exact commensur ability. In Fig. ^| we present 
new data for the activation energy for our longest sample, 
L = 1000a, with parameters W — 7a, Eq — 0.7K, and 
Ej = 0.9K. For each value of n, the activation energy Er 
was determined from the measured temperature depen- 
dence of the array resistance. The measurement was per- 
formed in the linear regime, at a small transport current. 
For this sample, the commensurate phase around the 
point no = l/3a exists in the domain |a| < olq ~ 0.009. 
The maximal value of the activation energy, Er « 12K, 
is reached at the commensurability point. Outside the 
Mott phase region, the resistance exhibits strong oscilla- 
tions; the activation energy vanishes almost periodically, 
with the period An = 1/L. We find two aspects of this 
data striking. 

Firstly, the regions of n corresponding to the Mott 
phase are extremely narrow (ac — 10~ 2 ). In the con- 
ventional picture, this would imply a weak interaction 
between the particles (compared to the one-particle band 
structure energies). Consequently, within the Mott phase 
the activation energies for particle transport must be also 
small. Quite contrary, the observed value of the resis- 
tance activation energy is about one order of magnitude 
larger than the energies Ec and Ej, which determine the 
single-vortex band spectrum. 
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FIG. 6. Activation energy of the resistance of an array con- 
sisting of 1000 x 7 cells with parameters Ej — 0.9K and 
Ec = 0.7K. The one-dimensional vortex densities n, no are 
measured in units of 1/a. The dashed line is a fit to the 
data, to extract the width of the Mott region. The cusp-like 
part of the figure corresponds to the Mott phase. Additional 
wiggles on that part may be related to boundary effects, see 
Section and Fig. |b. Inset: activation energy outside the 
Mott phase. The zeros of Ea(n) at n — no > 0.006 indicate 
the restored rigidity of the vortex chain. Note that in quali- 
tative agreement with Fig. the maxima of Er(ti) decrease 
with increasing \n — no|; the dotted line in the inset is a guide 
to the eye. 



(The resulting numerical value here refers to the parame- 
ters of the sample of Fig. ^). In the same approximation, 
the model form of the potential (^) becomes exact. With 
the help of Eq. (p3), the estimate (42) of the elastic con- 
stant can also be refined. For convenience, we give here 
the value of the product aK, instead of the value of the 
clastic constant: 



aK 



41n(V/a) 



3.5 x 10 4 K 



(54) 



The theoretical results ( |53|) and ( p4|) do not have any 
adjustable parameters, and are obtained within control- 
lable approximations. In contrast with this, Eq. ( p4j ) for 
the effective pinning cannot be used for the quantitative 
comparison with the experimental: The tight-binding ap- 
proximation we have used in Section III to estimate the 
suppression of the pinning potential is not applicable to 
the case Eq ~ Ej. (We note that Eq. (^3|) still allows one 
to reproduce the correct trend in the variation of ac with 
the ratio Ec/Ej). Therefore, we proceed in the following 
way. First, we find the soliton length x s from the experi- 
mental values of ac- Then, using the theoretical value of 
aK and the value of x s extracted from the data, we find 
the renormalized pinning potential Up S and the soliton 
activation energy E s . We will check that the renormal- 
ized pinning potential is indeed substantially lower than 
its bare value O.lEj. Finally, we will relate the found 
value of E s to the activation energy of the resistance for 
the experimental sample. 

The values of ac and a « 10~ 4 cm found experimen- 
tally allow us to estimate the soliton length from Eq. (147 



Secondly, the resistance R(n) exhibits strong oscilla- 
tions with the period An = 1/L outside the Mott region. 
These oscillations would not be expected in a model of 
almost-free quasiparticles within the delocalized phase. 

These two observations find a natural explanation in 
our model, which explicitly accounts for the long-range 
interaction forces between the vortices. 

In order to perform a detailed quantitative compari- 
son of the theory and experiment, first we improve the 
estimates ([H]) and ([l2|) for the interaction range A s and 
the elastic constant K respectively. In the experiment^, 
the contact bars were made by shorting the junctions 
at the lower and upper border of the array. We there- 
fore model the bars by superconducting strips of average 
width a « 10 _4 cm (the lattice constant of the array). 
The condition A s ^> W allows us to neglect non-local ef- 
fects in the solution of the magnetostatic problemEil, and 
to express A s in terms of the self- inductance L of the two- 
wire system. In addition, as the distance W between the 
wires exceeds significantly their width a, we can use the 
textbookB formula C = 41n(W/a). As a result, Eq. @ 
is replaced by 



167T 2 EjW\n(W/a) 



270a . 



(53) 



— = 890a w 0.09cm . 

ac 



(55) 



This length is really large. In fact, x s is about three times 
the length of the array L = 300a used in Ref . || to extract 
the activation energy En . This may explain why the val- 
ues of En found there are systematically lower than the 
activation energy for the longest array, see Fig. ^. A sin- 
gle soliton consists of about 300 vortices, and therefore 
its activation energy may exceed easily the single-vortex 
energy scales. Note also that x s exceeds considerably the 
interaction radius (|53| ) which gives us confidence in the 
applicability of the sine-Gordon equation ( [431). The ef- 
fective pinning potential, according to Eq. (|44|) , can be 
found as: 



2nx g 



— ~ 3.4 x 10" 3 K 
na 



(56) 



This energy is at least one order of magnitude smaller 
than its bare value 0.1-Ej, see Eq. (||). The reduction 
is apparently due to the quantum zero-point motion of 
individual vortices. Again, the strong suppression of the 
pinning energy guarantees the harmonic form of the pin- 
ning potential (|2^), and hence allows us to use the sine- 
Gordon equation for the solitons. Finally, using Eqs. (E3) 
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and pq), we find the soliton energy at the commensura- 
bility point: 



K 



(57) 



This energy exceeds significantly the boundary pinning 
energy. According to Eqs. (|l^) and (|50|), the latter con- 
tributes to Er less than 0.5K. We neglect this contri- 
bution, and therefore identify Er with the energy E s of 
the formation of a soliton. The calculated value (^) is 
somewhat lower than the measured Er. Still, we find 
the agreement quite impressive, having in mind the huge 
value of the elastic constant (|54|), calculated without any 
adjustable parameters. 

The large value of the elastic constant ([34]) results 
from the long-range nature of the inter- vortex interaction 
forces. In fact, the vortex chain in the incommensurate 
phase is so rigid, that Eq. ([52]) is inapplicable in the case 
of a sample only a thousand cells long. Away from the 
transition point (\a\ = cic), the elastic term in Eq. (52), 
which is supposed to be a small correction, is about 70 
times larger than the "main" boundary pinning term. 
Therefore, for the conditions of the experiment, the in- 
commensurate phase is well described by the model of 
a rigid vortex chain, see Section |0. This immediately 
explains the strong oscillations of the activation energy 
with the period An = 1/L in the incommensurate phase, 
see Eqs. @ and ©. 

There is a clear resemblance between the experimen- 
tally measured curve of Er(ji) (inset in Fig. ||), and the 
curve in Fig. ||, simulated with the help of Eqs. ( |l2| ) and 
(p6|). In agreement with the model of a rigid chain, the 
minima of Er{tl) reach zero at \n — no| ^ 0.005, and the 
maxima of Er(u) decrease with increasing |n— no\. The 
boundary term Eq. (12) has a maximum of sa 0.5K and 
accounts for the main contribution to the maxima of Er 
at n — no > 0.006. The vortex chain softens up only in a 
very narrow region around the transition point, so that 
the crossover region is of the order of An, see the inset 
in Fig. | 

To end this section, we reiterate that in the experi- 
ment the commensurability point no = 1 /3a was reached 
in the sample with W = 7a, which means the chain 
is stable against the formation of a zigzag structure at 
1/nW — 0.43. According to Eq. (|l5|), for a continu- 
ous system, the zigzag instability would already occur 
at 1/nW = 0.65, i.e., before the density n = l/3a is 
reached. Since the experimental data show no indica- 
tion of a qualitative difference between the arrays with 
W = 3a and W = 7a, we conclude that the array width 
W = 7a is narrow enough to allow suppression of the 
instability by the effects of discreteness. 



VII. DISCUSSION 

An external magnetic field applied to an array of 
Josephson junctions allows one to introduce vortices into 



it. A sufficiently weak field creates a linear chain of vor- 
tices in the quasi one-dimensional array. The ratio be- 
tween the periods of the vortex chain and the array of 
Josephson junctions is controlled by the value of the mag- 
netic field. The commensurate phase corresponds to the 
vortex analogue of a Mott insulator. Within this phase, 
the elementary excitation is a soliton consisting of a num- 
ber of individual vortices. The finite gap energy for the 
soliton translates into a finite activation energy of the re- 
sistance of the array. Each soliton transfers a fraction of 
the flux quantum through the array. In the incommensu- 
rate phase, the spontaneous proliferation of solitons and 
anti-solitons leads to the formation of a one-dimensional 
vortex liquid. This results in a finite vortex-flow resis- 
tance of the array. 

In this paper we have analyzed the commensurate- 
incommensurate transition for a one-dimensional vortex 
system in detail. The size and energy of the vortex soli- 
tons, which drive the transition, depend on two param- 
eters of the vortex chain. These parameters are: the 
elastic constant, and the pinning potential existing due 
to the discreteness of the array of Josephson junctions. 
The long-range nature of the vortex-vortex interaction 
leads to a large value of the elastic constant. On the 
other hand, the zero-point motion of each quantum vor- 
tex leads to a considerable suppression of the pinning 
potential. As a result, the size of the solitons turns out 
to be extremely large, about 300 vortices under the con- 
ditions of the experiments reported in this paper and in 
Ref. ||. This enables us to treat- the transition in the 
framework of the classical theoryEj. Our theory explains 
quantitatively the main experimental observations. 

We would like to conclude with the following remark: 
a quasi one-dimensional array of small superconducting 
islands connected by Josephson junctions can be used to 
study quantum phase transitions in two complementary 
ways. The first way relies on the control of the charge 
state of the islands by an external gate. In this case, 
a transition between the charge-localized and charge- 
delocalized phases can be observed in principle. The 
localized phase is a Mott insulator, with a finite gap 
for charge solitons, which play the role of elementary 
excitations. The delocalized phase behaves as a one- 
dimensional Luttinger liquid (see, e.g., Ref. ^ and refer- 
ences therein). The experimental observation of the two 
phases and the transition is difficult, as it is virtually im- 
possible to avoid the existence of random offset charges, 
which introduce strong disorder into the system. The 
other way is to study the commensurate-incommensurate 
transition in a system of vortices induced in the array by 
an external magnetic field (the case studied in this pa- 
per). This transition belongs to the same universality 
class as the Mott transition for charge derealization. A 
great advantage of the vortex system is that it is virtu- 
ally disorder-free. However, due to the large size of the 
solitons driving the transition, the critical region around 
the phase transition point is extremely narrow for the 
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arrays studied experimentally. To widen the critical re- 
gion, one should find a way to reduce the vortex-vortex 
interaction strength. That would open new possibili- 
ties of experimental investigations of the Luttinger liquid 
which is formed on the incommensurate side of the transi- 
tion. Trip-properties of the liquid are expected to depend 
cruciallyc.3 on the value of the fractional flux carried by 
the solitons. 
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